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Abstract. We use the Gutzwiller variational many-body theory to investigate 
the stability of orbitally ordered states in a two-band Hubbard-model without 
spin degrees of freedom. Our results differ significantly from earlier Hartree-Fock 
calculations for this model. The Hartree-Fock phase diagram displays a large 
variety of orbital orders. In contrast, in the Gutzwiller approach orbital order 
only appears for densities in a narrow region around half filling. 
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1. Introduction 

The investigation of orbital degrees of freedom has become an important field in 
theoretical solid-state physics over the past two decades. There are a number of 
materials which are believed to show phase transitions with orbital order parameters. 
Among them, the perovskite manganites, e.g., Lai_3.Sr3.MnO3, have attracted 
particular interest because of the colossal magnetoresistance behaviour which is 
observed in these materials. In theoretical studies on manganites, one often neglects 
the almost localised Mn t2 g orbitals and investigates solely the electronic properties 
of a systems with two e g orbitals per lattice site. In order to study the ferromagnetic 
phase of such a model, Takahashi and Shiba p] further neglected the spin-degrees of 
freedom because the Hund's rule coupling is assumed to align the spins in the two 
e g orbitals. 

The mean-field study in [1] found a surprisingly large number of stable orbitally 
ordered phases for the spinless e g model. However, it is well known that mean-field 
approximations tend to overestimate the stability of ordered phases in correlated 
electron systems. Therefore, the purpose of this work is to reinvestigate the two- 
orbital Hubbard model without spin-degrees of freedom by means of the Gutzwiller 
variational many-body theory. 

Our paper is organised as follows: The Hamiltonian and the different types of 
order parameters that we are going to investigate are discussed in section 2. In 
section 3 we introduce Gutzwiller wave functions and derive an approximate expression 
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for the variational ground-state energy. Numerical results are presented in section 4 
and a brief summary closes our presentation in section 5. 



2. Model system and types of order parameters 

We investigate a two-orbital (e g -type) Hubbard model [2] without spin-degrees of 
freedom, 

H =H + Hu, (1) 
X £ *ij'4*A, 6 , , (2) 

i.j b,b' — l 

H v = U^2 ni,in ii2 ■ (3) 

i 

b b' 

Here, the tight-binding parameters describe hopping processes between orbitals 
b, b' on cubic lattice sites R^ and R 3 , respectively. The Hamiltonian |T]) is formally 
equivalent to the standard one-band Hubbard model if the indices 6, b' are regarded as 
spins. However, the tight-binding parameters in @ would be unusual for a genuine one 
band model since they contain inter-orbital hopping-tcrms. The Hubbard parameter 
U in j3]) is derived from U = U 1 — J where V and J are the Coulomb and the exchange 
interaction between electrons in different e g orbitals. 

We restrict our investigation to systems with only nearest neighbour hopping since 
additional hopping terms would only destabilise the orbital order we are interested in. 
For e g orbitals \zi) = \x 2 — y 2 ), |z 2 ) = \3z 2 — r 2 ) and hopping parameters tdda = leV, 
tddS = OeV [3j, the one-particle Hamiltonian in momentum space reads 

Ha = 4;k,64;k,b' i (4) 



with [5] 



k 6,6' 



e l'M = ( cos ( k x) + cos (k y ) + 4 cos (k z ))/2 , 

2,2 



3(cos(fc :c )+cos(fc. y ))/2, (5) 
4ik = - v^cos (k x ) - cos (ky))/2 = e 2 £ . 

The k- integrated density of states D(e) that results from this band structure is shown 
in figure [TJ 

For our investigation of orbital order we need to introduce two more basis 
representations of the orbital space, 

W) = ^iW) + \Z2)), (6) 
\X2) = -^|(I Z 1> - l Z 2)), 



and 



\yx) = ^={\zx)+i\^)), (7) 
W) = — i| z 2> ) - 
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Figure 1. Density of states of the model system {4) 



:k • 



The dispersion relations in momentum space for the l x'- and the 'j/'-basis become 

4£ = ((2 - V3) cos (fc x ) + (2 + \/3) cos (fc„) + 2 cos (*»))/2 , (8) 
((2 + V3) cos (kx) + (2 - VS) cos (fc y ) + 2 cos (k z ))/2, 
cos (A; x ) + cos (fc y ) — 2 cos (k z ))/2 = 
cos (fc x ) + cos (fc. y ) + cos (k x ) = , 
((1 + V3i) cos (k x ) + (1 - V3i) cos (k y ) - 2 cos (k z ))/2 , 

The interaction term ([3|) has the same form for all three basis representations 
£ = x,y,z. 

We describe orbital order in our model system {!]) through the parameters 



2.2 
E ic;k 



1.1 

1,2 

2,1 



(9) 



C ltb C i,Zb ' 



,»/2, 



(10) 
(11) 



for each of the three representations £. By using the Pauli matrices 1, f x , f y , f z the 
order parameter can also be written as (compare reference [1]) 

1 

2 



6,6' 



(cL h (^)6,6'C ?; . ) 



(12) 



Besides the orbital character of the order-parameter we need to specify its lattice 
site dependence. Following reference [T| we consider orders of the form 

t^i = T£,o exp (iQRj) (13) 

with commensurate vectors Q which belong to the T point (Q = (0, 0, 0)), the R point 
(Q = (7r,7r,7r)), the X point (Q = (0,0, n)), and the M point (Q = [k,tt,0)). The 
real parameter r^ o in (p~3|) is independent of the lattice site vector and assumed 
to be positive. Note, that for vectors Q ^ (0,0,0) equation (|13|) divides the lattice 
into an '.A '-lattice with a majority £i occupation (r^ ; , > 0), and a 'B'-lattice with a 
majority £2 occupation (t£ ; , < 0). 
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3. Gutzwiller wave-functions 

3.1. Definitions 

For an investigation of the Hamiltonian ((T|) we use Gutzwiller variational wave 
functions [4] which are defined as 

i* g ) = IP<i*°> • ( 14 ) 

i 

Here, |Wo) is a normalised quasi-particle vacuum and the local correlator has the form 

Pi - J~] h,irhi,i (15) 
i 

where rhj = |/)(/| projects onto the four local configuration states |J), i.e., the empty 
state |0), the doubly occupied state \d) and the two single electron states |£i),|£2)- 
Note, that and |^) depend on the orbital representation £ whereas the states |0) 
and \d) are invariant under the orbital transformations ([6]) and ([7]). For the variational 
parameters Xij we make an Ansatz which is consistent with the spatial symmetry of 



the order parameter, 

^i,0 = ^0 , (16) 

Xi,d = X d , (17) 

A i>6 =X S + 5X S exp (iQRj) , (18) 

Xifo =X S - SX S exp (iQRi) , (19) 



where the parameters Ag,Arf, X S ,6X S are independent of the lattice site vector R;. 

A more general ansatz for the Gutzwiller correlator P which also includes non- 
diagonal configuration transfer operators |/)(/'| has been studied in reference [5j. As 
we will see below (section 13. 2p , it is sufficient for our calculations in this work to 
consider only diagonal operators in Pj. 

3.2. Evaluation in infinite dimensions 

For any practical use of a variational wave-function it is essential that the expectation 
value of the Hamilton can be calculated. However, despite the simplicity of the 
Gutzwiller wave-function the evaluation of 

£var = (H)v a (20) 

poses a difficult many-particle problem that cannot be solved in general. Gutzwiller 
introduced an approximate evaluation scheme that was based on quasi-classical 
counting arguments. More recent derivations of this approximation can be found 
in references [51 E]. Analytically exact evaluations were later found to be possible 
in one dimension [8l [9l I10j . and in infinite spatial dimension [llj . The results of 
the latter evaluation turned out to be equivalent to the Gutzwiller approximation 
for systems which can be studied within this approach. The evaluation scheme in 
infinite dimensions was later generalised for the investigation of multi-band Hubbard 
models [H] and superconducting systems [SJ Q3] . We will use these exact results in 
infinite dimensions as an approximation in order to evaluate expectation values of our 
Hamiltonian |T]). 

We only consider single-particle wave functions |^o) i n 03) f° r which the local 
density matrix C° with the matrix elements 

^,^(^/«>. E 4v (2i) 
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is diagonal with respect to b, b' . Finite non-diagonal elements in C° could only appear 
if we were mixing different order parameters. 

As shown in [5] the four parameters A, / for a lattice site i have to obey the 
constraints 

1 = <A 2 )*o , (22) 

CSfc.e, = <W,fc W*o ■ ( 23 ) 
Our correlation operator ([151) automatically fulfils the constraints ([23| for b ^ b'. This 
is the reason why it was allowed in the first place to include only diagonal operators 
riiij in p5]) . Consequently instead of ([25)) we only need to consider the diagonal 
constraints 



o _ / D 2 



= <*?ftiA>*o • ( 24 ) 
All constraints can be solved explicitly if we use the results for local expectation values 

rriij = (rhij)y G = \ 2 lI {m iJ )q lQ (25) 

which hold in infinite spatial dimensions. With (|25[) we can use the expectation values 
rriij as new variational parameters instead of The constraints then read 

1 = m,j + m 4 , d + m, A + m,,f 2 , (26) 

n it b = TO ^«f + TO M ■ (27) 
and can be readily solved by expressing all local occupancies in terms of the average 
numbers of doubly occupied sites m^d, 

m i; = 1 - n° ?1 - + ra M , (28) 

m *,«b = "i.ft, ~ m M • (29) 
Apart from the still unspecified one particle wave function j^o) and the corresponding 
local density matrices Cf the probabilities m^d are the only remaining variational 
parameters. Note, that leads to 

= («t,&)*G = ( 30 ) 

for the orbital densities in the correlated Gutzwiller wave-function. We skip the 
explicit declaration of the orbital representation £ for the rest of this chapter and 
write b instead of m &U indices. 

In infinite dimensions the expectation value of the one particle Hamiltonian @ 
becomes [T2] 

E = (#o)* G = Yl ]C Vli,bqj,b'ti,j (c!,&Cj,6')*o = (#o)*o (31) 
i.j b.b' 

where we introduced the well known Gutzwiller loss factors [4] 

1 2 

= o/i OA (y/ m i,Q m i,b + ( 32 ) 

and b is defined via 1 = 2 and 2 = 1. 

The lattice symmetry of the order parameter leads to further simplifications. We 
introduce the majority and the minority orbital density 

n°± = n° ± t (33) 

such that 

<i = n\ , < 2 = n% . (34) 
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The upper and lower signs in (|34[) . and in corresponding equations below, belong to 
lattice sites i £ A and i £ B, respectively. For the other local expectation values we 
find 

mi,d = m d , (35) 



m, 
m, 
m. 



1 = m± = n° ± - m d , (36) 

2 = m T = - m d , (37) 
= 1 — 2n° + m d . (38) 

A similar notation is introduced for the g-factors 

Qi,i = Q± , ?i,2 = Qt ( 39 ) 

where 

I 2 

q± = ^oTl ^oT (V m TO ± + V m T m d) ■ (40) 

n ± [L — n ± ) 

The expectation value (|31[) splits into four components 

E =q + E ++ +q_E— +^qI(E+- + E~+) (41) 

which belong to the four different hopping channels between majority ('+') and 
minority ('— ') states. In momentum space, the one-particle expectation values in 
l} can be written as 

E " J = lEE**' (4((i^>*. +^6 > 6'(4+Q >6 4,6'>*c (42) 
k 6,6' 

+ w ' M 6 3 6'(4,6 6 k+Q,6')*o + ^'^4 4 6'(c k +Q,6 6 k+Q,6')*o) 

"i 

elements of the matrices 



where u, u>' represent the + or — signs and the coefficients M7 b , are given as the 



I I) , M^(\ \ \, (43) 



M = U -i J < M - v-i i 

For the evaluation of the expectation values in (|4"2j) we need to determine the one- 
particle wave function |^o)- Following reference [SJH3], |\I>o) is given as the ground 
state of the effective one-particle Hamiltonian 

Hf = ff - r^exp (iQR0(«,,i - ni, 2 ) , (44) 

i 

where Hq was introduced in (|3ip and the term proportional to the variational 
parameter r\ allows to vary the order parameter to- 

In this work, we only aim to investigate the stability of the orbitally unordered 
state. Therefore we just need to analyse the energy expression (j41j) for small values of 
the order parameter To. An expansion of the g-factors f)40[) up to second order in To, 

q± w q Q ± qxTo + q~2T% (45) 

yields 
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Note, that the coefficients q 7 are still functions of n° and m d . To leading order in t 
the effective Hamiltonian (|4"4")) becomes 

Hf = VoJ2J2 e k b 4,&Ck,6' - V J2 (4, A+q,i - 4, 2 c k +Q,2) (47) 

k 6,6' k 

since we can set g^6 = <?o i n H' Q , cc l- (EU- The one-particle Hamiltonian (|47|) is easily 
diagonalised numerically. This diagonalisation yields the coefficients in the quadratic 
expansion of (|42[) 

£ ±± = E + E 2 tI , (48) 
E±t = E ' + E' 2 t* , (49) 
and, consequently, of the variational ground state energy 

E vax = q Q El ot + (qoE^ + q 2 E^ - E' ^j t 2 + Um d . (50) 

Here, we introduced 

E^ = 2(E 1 + E' y ) (51) 

for 7 — or 2. The minimisation of £„ ar with respect to m,i can be carried out for 
ro =0 



E 



tot 9 «0 



dm d 



+ U = (52) 

and determines the optimum number rhd of doubly occupied sites as a function of U 
and n°. This allows the expansion of the variational energy purely in terms of tq 

£var = 90^ Ot + C ro T 2 . (53) 




Figure 2. Critical interaction strength in Hartree Fock-theory for the orders 

(straight), rf, (dashed), Tp (long-dashed), t?j (dotted), and rj^ (dash-dotted). 
These phases are not stable within the Gutzwiller theory. 
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A negative sign of the coefficient 

c T0 = q El ot + q 2 El ot -E'S= Or, (U, n°) (54) 

in (|53)l indicates the instability of the unordered state. Note that a positive c Tfl not 
necessarily proves the stability of the unordered state since it does not exclude first 
order transitions. 

In Hartree-Fock theory a quadratic expansion of the ground-state energy leads to 
E™. = E«* + {E™ - C/)r 2 , (55) 
and the critical interaction strength in this approach is therefore given as 

U$ F = E%* . (56) 



4. Results 

In figures [2] and [3] we show the critical interaction strength (|56|) in Hartree-Fock theory 
as a function of density for the various types of orbital order introduced in section [2] 
Our data agree very well with those reported in reference pQ. Note that in Hartree- 
Fock theory the critical interaction strength is finite for all densities n° > and 
diverges only in the limit n° — > 0. 

The phases in figure [2] are not stable within our correlated Gutzwiller approach 
for all densities and interaction parameters. This holds in particular for the order 
parameter which has surprisingly small critical values U^ F around quarter filling 
2n° « 0.5. 

For the four surviving phases (figure [3]) we show the ratio Uq W /Uq F of the 
critical parameters in Gutzwiller and Hartree-Fock theory as a function of density in 




Figure 3. Critical interaction strength in Hartree Fock-theory for the orders 
T r(X) (straight), (dashed), (dotted), and -rj (dash-dotted). 




Figure 4. Critical interaction strength in Gutzwillcr-thcory for the orders t±( , 

1 ( A ; 

(straight), t^ x ^ (dashed), (dotted), and rj^ (dash-dotted). The arrows 

indicate the critical densities 2n[? for U — * oo. 



figure 2J Apparently, there is only a narrow window of densities around half filling 
2n° ~ 1 where orbital order occurs in the Gutzwiller theory This is in stark contrast 
to the Hartree-Fock findings in figure [3] For three of the phases the orbital order may 
disappear again if U exceeds some second critical value Uq w > Uq w . This behaviour 
is also different from the Hartree-Fock theory in which the orbital order is stable for 
all U > Uq F . Mathematically, the appearance of the second critical parameter U^ w 
is due to the fact that qi in (|45[) has a minimum as a function of mj. Hence, the 
coefficient c To may have one or two roots as a function of U , depending on the other 
parameters in (|54|) . The appearance of the second transition is therefore a genuine 
many-particle effect. Only one of the four orders, rj^ , is unstable in the limit U — ► oo 
for all densities. For the three other orders we find critical densities n° with a stable 
order for all n° > n° in the limit U — > oo. The values of n° are displayed by arrows 
in figure [4j 

The failure of Hartree-Fock theory to describe the orbitally ordered phases of 
the two-band model is not surprising. It is well known, for example, that Hartree- 
Fock theory also grossly overestimates the stability of ferromagnetic ground states 
in the one-band Hubbard model |15j . Ferromagnetism in this model requires very 
peculiar densities of states or very large local Coulomb interactions |16j . The simple 
Stoner criterium of the Hartree-Fock theory, however, predicts ferromagnetism for all 
densities and for small Coulomb interactions. 

Orbital order in the two-band model at infinite U has also been investigated in 
reference [17] . In their work the authors report a complete disappearance of all types 
of orbital order seen in the Hartree-Fock theory. 
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5. Summary 

We have reported results for orbital order in a two-orbital Hubbard-modcl without 
spin-degrees of freedom. In a previous work which was based on Hartrce-Fock 
calculations, this system seemed to exhibit a surprisingly large variety of orbitally- 
ordered phases. In our study in this work we have employed the Gutzwiller variational 
theory which is known to be more reliable than the Hartree-Fock theory for systems 
with medium to strong local Coulomb interaction. Most of the phases found in the 
Hartree-Fock approach turned out to be unstable in the Gutzwillcr-thcory. Orbital 
order only appears for densities near half filling in our calculation. Unlike in Hartree- 
Fock theory, it may happen that an orbitally ordered phase which is stable for 
correlation parameters U > Uq W becomes unstable again if U exceeds a second critical 
value Uq W > Uq w . This second transition is a genuine many-particle effect. 

Our findings show that the stability of phases with broken symmetry for correlated 
electron systems can be grossly overestimated by the Hartree-Fock mean-field theory. 
It is quite likely that similar problems appear in LDA+U calculations where the local 
Coulomb interaction is also treated on a Hartree-Fock level. 
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